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Abstract 



In this paper, we study a linearized two-dimensional Euler equation. This equa- 
tion decouples into infinitely many invariant subsystems. Each invariant subsystem is 
shown to be a linear Hamiltonian system of infinite dimensions. Another important 
invariant besides the Hamiltonian for each invariant subsystem is found, and is utilized 
to prove an "unstable disk theorem" through a simple Energy-Casimir argument [1]. 
The eigenvalues of the linear Hamiltonian system are of four types: real pairs (c, — c), 
purely imaginary pairs (id, —id), quadruples (±c ± id), and zero eigenvalues. The 
eigenvalues are computed through continued fractions. The spectral equation for each 
invariant subsystem is a Poincare-type difference equation, i.e. it can be represented as 
the spectral equation of an infinite matrix operator, and the infinite matrix operator is 
a sum of a constant-coefficient infinite matrix operator and a compact infinite matrix 
operator. We have obtained a complete spectral theory. 
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I. Introduction 

In [2], Henshaw and Kreiss numerically studied the propagation of perturbations in the 
solution of the two dimensional incompressible Navier-Stokes equations with no body force, 
and at high Reynolds numbers. They numerically solved the equations starting with some 
smooth initial data for which the Fourier modes have random phases. This flow evolved 
over time, initially through a complicated state containing many shear layers and then into 
a state of large vortex structures (with shear layers) which persists for long time. 

1. Changing the viscosity, the large vortex structures do not change much. 

2. Adding high mode perturbation to the initial data, the large vortex structures do not 
change much. 

3. Changing the Laplacian A to A 2 , the large vortex structures do not change much. 

4. Changing the Laplacian A to an operator which is A for low modes, and is A 2 for high 
modes, the large vortex structures do not change much. 

In [3], Matthaeus et al. studied the same problem, and found similar results. Matthaeus et 
al. run the numerics for much longer time. These numerics indicate that for relatively long 
time (not infinite long time), the solution to the 2D N-S equation (without body force, i.e. 
decaying turbulence) has the large vortex structures. Such structures persist in the solution 
to the 2D Euler equation by the claims 1, 3, 4 above. (Cf: Under decay boundary conditions, 
the Kato's theorem states that for finite time the solution to 2D N-S equation converges to 
the solution to 2D Euler equation in norm as viscosity approaches zero, see [4].) Moreover, 
the large vortex structures are stable with respect to the change in initial data. 
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In [5], Robert and Sommeria studied the organized structures in two-dimensional Euler 
fluid flows by a theory of equilibrium statistical mechanics. The theory takes into account 
all the known constants of motion for the two-dimensional Euler equations. The microscopic 
states are all the possible vorticity fields, while a macroscopic state is defined as a probability 
distribution of vorticity at each point of the domain, which descibes in a statistical sense 
the fine-scale vorticity fluctuations. The organized structure appears as a state of maximal 
entropy, with the constraints of all the constants of motion. The vorticity field obtained as 
the local average of this optimal macrostate is a steady solution of the Euler equations. 

The above numerical results show that certain relatively long time large vortex struc- 
tures for 2D Euler equation persist for 2D N-S equation at high Reynolds numbers. Such 
structures are stable with respect to the change in initial data. We believe such structures 
are the exhibition of certain unstable manifolds. The above theoretical results show that 
the probabiltity mean of such structures for 2D Euler flows are steady solutions to the 2D 
Euler equations. Thus, we believe that certain unstable manifolds, if there is any, of steady 
solutions to 2D Euler equations are responsible for such large vortex structures. Therefore, a 
dynamical system study on certain unstable manifolds of certain steady solutions to 2D Euler 
equations, and their persistence for 2D N-S equations at high Reynolds numbers is crucial 
for studying the large vortex structures. Indeed, we have built a chaos-molecules-model on 
2D turbulence based upon the above motivations [6]. We believe that our chaos- molecules- 
model captures the qualititive frames of the hyperbolic structures, and the energy inverse 
cascade and the enstrophy cascade nature for 2D turbulence. 

In this paper, we study the linearized 2D Euler equation at a fixed point. This study is 
the base for future analytical studies on the unstable manifolds for the 2D Euler equation. 
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In [7], we have begun numerical studies on the unstable manifolds for the 2D Euler equation. 

The current study is also important in the linear hydrodynamic stability theory. By 
utilizing Energy- Casimir method [1], we obtain an unstable disk theorem which is not in the 
category of the classical Rayleigh theorem [8] . 

Next we discuss the approaches used in this study. Through the Energy- Casimir method, 
nonlinear stabilities of various types of two-dimensional ideal fluid flows have been established 
[9] [10] [11] [1]. Below we give a brief description on the Energy- Casimir method. Let D be 
a region on the (x, y)-plane bounded by the curves Tj (i = 1,2), an ideal fluid flow in D is 
governed by the 2D Euler equation written in the stream-function form: 

J^=[W,VA^] , (1.1) 

where 

ryw, vA^l = ^±^t _ ; 

dx dy dy dx 

with the boundary conditions, 

. d r dip , 

For every function f(z), the functional 

F = J j D f(A^)dxdy (1.2) 
is a constant of motion (a Casimir) for (1.1). The conditional extremum of the kinetic energy 

E = \l / D V ^- V ^ dxdy ( L3 ) 
for fixed F is given by the Lagrange's formula [9], 

5H = 8(E + XF) = , = Xf (Ai/j ) . (1.4) 
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where A is the Lagrange multiplier. Thus, i/j is the stream function of a stationary flow, 
which satisfies 

ij = $(A^ ) , (1.5) 
where $ = \ f. The second variation is given by [9], 

5 2 H = \j J d { V0 • V0 + $'(A^ ) (A0) 2 Jcfecfy . (1.6) 

Let ip = ipo + V 9 be a solution to the 2D Euler equation (1.1), Arnold proved the estimates 
[10]: (a), when c < &(Aip ) <C,0<c<C <oo, 

J jv^(f) • V<p(t) + c[Ay?(t)] 2 } cfccfy < J J {v<^(0) • V<^(0) + C[Av?(0)] 2 j cfecfy, 

for all t e (-oo, +oo), (b). when c < -$'(A^o) <C, 0<c<C<oo, 

J J {c[Ay?(t)] 2 - Vy?(t) • Vy?(t)j dxdy < / J |c[Av?(0)] 2 - V^(0) • V^(0) J dxdy, 

for all t G (— oo,+oo). Therefore, when the second variation (1.6) is positive definite, or 
when 

J J d |V0 • V0 + [max$'(A^ )] (A0) 2 | dxdy 

is negative definite, the stationary flow (1.5) is nonlinearly stable (Liapunov stable). In this 
paper, we have found an invariant for the linearized 2D Euler equation, and use this invariant 
together with an Energy-Casimir type argument to study linear stability, and to prove an 
unstable disk theorem. The linearized 2D Euler equation is an infinite dimensional linear 
Hamiltonian system. For finite dimensional linear Hamiltonian systems, it is well-known 
that the eigenvalues are of four types: real pairs (c, — c), purely imaginary pairs {id, —id), 
quadruples (±c±id), and zero eigenvalues [12] [13] [14]. The same is true for the linearized 
2D Euler equation. The eigenvalues are computed through continued fractions following the 



work of Meshalkin and Sinai [15]. The linearized 2D Euler equation can also be written 
in an infinite matrix form. The spectral equation of the infinite matrix operator defines 
a Poincare-type difference equation [16] [17]. That is, the infinite matrix operator can be 
written as the sum of a constant-coefficient infinite matrix operator and a compact infinite 
matrix operator. In this paper, we follow a spectral theory developed by Duren [18] to 
study the spectra of the constant-coefficient infinite matrix operator through characteristic 
polynomials. Then we apply the Weyl's essential spectra theorem to the perturbation of the 
constant-coefficient infinite matrix operator by the compact infinite matrix operator [19], to 
achieve a complete spectral theory. 

Finally, we discuss some preliminaries on the 2D Euler equation. Consider the two- 
dimensional incompressible Euler equation written in vorticity form, 



dt dx dy ' 



du dv 
dx dy 



(1.7) 



under periodic boundary conditions in both x and y directions with period 2n, where Q is 
vorticity, u and v are respectively velocity components along x and y directions. We also 
require that both u and v have means zero, 



I I u dxdy = v dxdy = 0. 

Jo Jo Jo Jo 

Expand Q into Fourier series, 

fcgZ 2 /{0} 

where u^k — , k — (/ci,/c 2 ) T , X = (x,y) T . In this paper, we confuse with (0, 0) T , the 
context will always make it clear. By the relation between vorticity Q and stream function 



dv du 



dx dy 

where the stream function \1> is defined by, 



d^ _ d^ 

dy ' dx 



the system (1.7) can be rewritten as the following kinetic system, 



k=p+q 



where A(p, q) is given by, 



A(p,q) = ^[\q\ 2 -\p\ 2 }(piq2 ~ P2qi) 



\q\- 2 -\p[ 



where \q\ — q\ + q\ for q — (qi, q 2 ) , similarly for p. 



Pi Qi 

P2 <?2 



(1.8) 



(1.9) 



Remark 1.1 Notice that direct calculation shows that the nonlinear term in (1.8) is J2k= P +q A(p, q) uj p uj q ^ 
where A(p,q) = \p\~ 2 (qiP2 -Piq.2); then, A(q,p) = \q\~ 2 (piq2 -QiP2)- The A(p,q) in (1.8) is 
the average ofA(p,q) and A(q,p), A(p,q) = \[A(p,q) + A(q,p)]. 



For any two functionals F 1 and F 2 of {uo k }, define their Lie-Poisson bracket: 



k+p+q=0 



Qi Pi 

<?2 P2 



dF\ 0F2 
duJp~ duJg' 



(1.10) 



Then the 2D Euler equation (1.8) is a Hamiltonian system [20], 



= {^k, H}, 



(1.11) 
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where the Hamiltonian H is the kinetic energy, 

H = \ £ |fc|-> fe | 2 . (1-12) 
z feezV{o} 

Following are Casimirs (i.e. invariants that Poisson commute with any functional) of the 
Hamiltonian system (1. 11): 

J n = u kl ---u kn . (1.13) 

fci+-+fc„=0 

The following Proposition is concerned with the equilibrium manifolds of the 2D Euler 
equation (1.8). 

Proposition 1 For any k G Z 2 /{0}, the infinite dimensional space 

El = {{ou k ,} u k , = 0, if k' + rk, Vr G , 
and the finite dimensional space 

E 2 k = {{cu fc ,} I cu k , = 0, if\k'\ , 
entirely consist of fixed points of the system (1.8). 

Proof: Let A; G Z 2 /{0}, {^} G For any p, q G Z 2 /{0}, ^ implies that p = 

rik°, q = r 2 k° for some r±, r 2 G -R; then, A(p, g) = 0; thus we always have A(p, q) oOpOO° q = 0, 
and {cu° k } is a fixed point of (1.8). Let k° G Z 2 /{0}, {u° k } G For any p, q G Z 2 /{0}, 

cjpCj° 7^ implies that \p\ = \k°\, \q\ = \k°\; then, A(p,q) = 0; thus we always have 
A(p, q) u)pU® = 0, and {w k } is a fixed point of (1.8). □ 

Fig.l shows an example on the locations of the modes (k' = rk) and (\k'\ = \k\) in the 
definitions of E\ and E\ (Proposition 1). 

The paper is organized as follows: Section II is on the formulations of the problem. 
Section III is on the Liapunov stability. Section IV is on the properties of the eigenvalues 



of the linearized 2D Euler equation as a linear Hamiltonian system. Section V is on the 
continued fraction study of the eigenvalues. Section VI is on the infinite-matrix study of the 
spectra of the linearized 2D Euler equation. Section VII is the conclusion. 
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II. The Formulations of the Problem 

Denote {uJk}kez 2 /{o} by u. Consider the simple fixed point u*: 

lo* p = T, u* k = 0, if A; or -p, (II.l) 

of the 2D Euler equation (1.8), which belongs to the two-dimensional intersection space 
Ep fl Ep (Proposition 1), where T is an arbitrary complex constant. The linearized two- 
dimensional Euler equation at u* is given by, 

io k = A(p, k - p) T co k ^p + A(-p, k + p) f cuk+p ■ (II.2) 

This is the linearized two-dimensional Euler equation that we are going to study in this 
paper. 

Definition 1 (Classes) For any k £ Z 2 /{0}, we define the class to be the subset of 
Z 2 /{0}: 

S^. = | A; + np G Z 2 /{0} n G Z, p is specified in (II. 1)^. 



See Fig. 2 for an illustration of the classes. According to the classification defined in Definition 
1, the linearized two-dimensional Euler equation (II. 2) decouples into infinite many invariant 
subsystems: 

^k+np = A(p,k + (n-l)p)T u; k+{n _ 1)p 

(II.3) 

+ A(-p, k + (n+ l)p) T u k+{n+1)p . 

Each invariant subsystem can be rewritten as a linear Hamiltonian system as shown below. 
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Definition 2 (The Quadratic Hamiltonian) The quadratic Hamiltonian H.- k is defined 



as: 



Hi 



-2 Im{^> n r A(p,k + (n-l)p) 

n<=Z 



U k+(n-l)p U k+np 



(II.4) 



Pi h 
P2 fa 



U k+(n-l)p ^k+np 



\m\ r Pn Pn-1 

where p n = [\k + np\~ 2 — \p\~ 2 }, " Im " denotes " imaginary part ". 



Then the invariant subsystem (II. 3) can be rewritten as a linear Hamiltonian system, 

1 U k+np ~ Pn q- 



(II.5) 



k+np 



Let 



^k+np = a n + i(3n, ne Z , 



i.e. a n and (5 n are the real and imaginary parts of u) k+np . Then the linear Hamiltonian 
system (II. 5) can be rewritten in the form, 



a. 



2 Pn ap n ' 



(II.6) 



8 = _I o-i 



where 



fa Pi 
fa P2 



Pn Pn-1 

neZ 



Y r (oi, n (3 n —\ o; n _i/9 n ) 
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where r = T r + i Tj, i.e. T r and Ij are the real and imaginary parts of T. Explicitly, 

Pn+l Tr Oin+l ~ Pn-1 Tr OL n -\ 



a. 



h Pi 
h Pi 



Pn 



Pn+l Tj Pn+l + Pn-1 Tj 

Pn-1 Tr Ai-1 — Pn+l IV Az+1 



(II.7) 



fci pi 
h P2 



+Pn+1 Tj « n+ i + Pn-1 r? «n-l 

If we rescale the variables as follows, 



1 



T 



fa Pi 
fa P2 



t , (X n p n Oi n , P n p n P n , 



The linear Hamiltonian system (II. 7) can be rewritten in the simpler form, 



da n /dr 



dp n /dT 



r r (a n +i - a„_i) + r, (p n+ i + 



-Pn 



T r (P n -i - Pn+i) + Ti (a n +i + a n -i) . 
Next we discuss the fixed point u* from a variational-principle point of view. Consider the 
kinetic energy and the enstrophy of the 2D Euler equation (1.8), 

E = \ \k\~ 2 \uj k \ 2 , J= \ u k\ 2 ■ 

If we extremize the kinetic energy E for fixed enstrophy J = c (a constant), we have the 



critical states by the Lagrange formula, 

J - c = , 



dL 



dL 

du> k 



2 + \)cu k = 0, VA; e Z 2 / {0} , 
where L = 2E + X(J — c). Thus, the critical states satisfy the relations, 

' A = - \q\~ 2 , for some q G Z 2 /{0}, 

u k = 0, if \k\ ^ \q\ , 



\u k \ =c. 

\k\ = \q\ 



(11.8) 
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Thus, we have the critical manifold for extremizing the kinetic energy for fixed enstrophy, 



= lu u k = 0, if \k\^\q\, K| 2 = c f> 

L \k\=\q\ J 



which is a sub manifold of the equilibrium manifold E^ (Proposition 1). Denote by / the 
linear combination, 



I = 2E- \p\~ 2 J. 



(II.9) 



Then, we have the following variational principle for the fixed point oo* (II. 1). 

Variational Principle: The fixed point oo* is a conditionally critical state of the kinetic 
energy E for fixed enstrophy J = |T| 2 , and is an absolute critical state of I. 

III. Liapunov Stability 

Definition 3 (An Important Functional) For each invariant subsystem (II. 3), we de- 
fine the functional I k which is the restriction of the functional I (II. 9) to the class 



*k -^(restricted to s^) 



= E^+^r 2 -br 2 } 



(III 1) 



^k+np 



Lemma III.l I k is a constant of motion for the system (II. 5). 



Proof: Differentiating Iy we have 



U k+np U k+np + ^k+np^k+np 



-i E 

n&Z 



k+np ' ~~~ U k+np a 

UUJ k+np UUJ k+np 
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Pi h 
Pi h 



E U k+np\ V PnPn~l^k+(n~l)p 



-Fp n Pn-lVk+(n- 



l)p 



0. 



This completes the proof of the lemma. □ 

Next we define the concept of disk in Z 2 / {0} which is needed in the unstable disk theorem 
to be proved below. 

Definition 4 (The Disk) The disk of radius \p\ in Z 2 / {0}, denoted by D\ p \, is defined as 



D lpl = [keZ 2 /{0} 



\k\ < \p\ 



The closure of D\ p \, denoted by D\ p \, is defined as 

D\ p \ = \kez 2 /{U} | \k\<\ P \ 

See Fig. 2 for an illustration. Next we prove the unstable disk theorem using a simple Energy- 
Casimir type argument [20] [1]. 

Theorem III.l (Unstable Disk Theorem) // Hj.C\D\ p \ = 0, then the invariant subsys- 
tem (II .3) is Liapunov stable for all t e R, in fact, 



E ^k+npit) < ° E K+np(°) 
neZ neZ 



Vt e R, 



where 



max{-p n } 



mm{-p n } 

n&Z 



-1 



, < a < oo 
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Proof: By lemma III. 1, 1^ is a constant of motion for the invariant subsystem (II. 3); then 



2 2 

E Pn U k+np (t) =Y,Pn ^fe+n P (°) > Wt E R . 
neZ neZ 



(III.2) 



If n D\ p \ = 0, then 



\k + np\ > \p\ , Wn <E Z , 



Thus, there exists a constant 5 > 0, such that 



<5 < -p n < 2 . 



(III.3) 



By (III.2), 



min{-p n } E k +np (0 



< 



E Pn U k+np (t) 



n&Z 



~T,Pn ^+np(°) 



neZ 



< max{-p n }^L +?!p (0) 



neZ 



n&Z 



that is, 



E K+npW < ° E | W fc+np(°) 
n6Z neZ 



where 



max{-p n } min{-p„} 



-1 



By relation (III. 3), 



-5<a< 25- 1 . 
2 



This completes the proof of the theorem. □ 



16 



Remark III.l // k \\ p, i.e. 3 real scalar a such that k = a p , then the invariant subsystem 
(II. 3) reduces to 

thus, it is obviously Liapunov stable for all t G R, in fact, this is a linearization inside the 
equilibrium space E p of the 2D Euler equation (cf: Definition 1). 

If \k\ — then the invariant subsystem (II. 3) decomposes into two decoupled systems, 
^k+np = A (p, k + (n- l)p) Tu k+{n _ 1)p 



+A(-p, k + (n + l)p) TuJ k+{n+1)p , (n > 1), 



(III.4) 



CO 



k+np = A (P, k + (n- l)p) ToJ U(n _ 



+A(-p, k + (n + l)p) TuJ k+(n+l)p , (n < -1), 
where A(p, k) = 0. The equation for uo k is 

u k = A(p, k-p) Yu k _ p + A(-p, k + p) Tu Up . 
Each of (HI. 4) and (III. 5) is a Hamiltonian system with the Hamiltonian, 



(III.5) 



(III.6) 



Pl 


k 


Im | 


P2 


k 





Im I Yl V PnPn~lU" k+{n _ l)p UJ~ k+np , 



.71=1 



Pi k 



Im I J2 V PnPn-iu k+{n _ iy! p k+np \ ; 



n=-l 
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which has the same representation as (II. 5). Denote by l£ and 7- the restrictions of If, to 
the systems (III.4) and (III.5): 



^ - E Pn l^k+np 



n=l 



h = E Pn 

n=-l 



U k+np 



(III.7) 
(III.8) 



Lemma III. 2 If \k\ = \p\, then l£ and 7- are respectively constants of motion for the 
systems (III. 4) and (III. 5). 



Proof: The same with that for Lemma III. 1 . □ 

Theorem III. 2 (Half Class Stability Theorem) If \k\ = \p\, k + p D\ p \, then the 
linear Hamiltonian system (III. 4) is Liapunov stable for all t G R, in fact, 



Ek+np(*) <*Ek+np(°) . v * e ^ 



n=l 



where 



max{-p n } 



n>l 



n=l 



min {-p n } 

n>l 



, < er < oo . 



If \k\ = \p\, k — p ^ D\ p \, then the linear Hamiltonian system (III. 5) is Liapunov stable for 
all t G R, in fact, 



E k+™(*) ^ E h +n M , vtGi?, 



n=-l 



n=-l 



where 



max{-p n } 

n< — 1 



i -1 



min {-p n } 

n< — 1 



< a < oo . 



Proof: The same argument on 1^ and 7- as that on 7^ in the proof of Theorem III. 1 . □ 
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Remark III. 2 If k — (p2,—Pi) T or k = (—p2,Pi) T , then k + p D\ p \ and I z — p D\ p \. 
Therefore, by Theorem III. 2, both systems (III. 4) an d (III. 5) are Liapunov stable for all 
t G R. Points in are on the tangent lines to the circle of radius \p\ at k (cf: Fig. 2). 
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IV. Properties of the Point Spectrum for the Linearized| 
Two- Dimensional Euler Equation as a Linear Hamil-| 
tonian System 



In this section, we study the properties of the eigenvalues for the linear Hamiltonian system 
(II. 6). The right hand side of (II. 6) defines a linear operator denoted by C, i.e. 



a 

P 



V 



(IV.l) 



where 



a = 


(• 


• • a. 


-l «o £*i 




P = 


{•••P-iPoPi- 


••) T ; 


t = 


(• 






•) T , 


V = 


{■■■V-i Vo Vi- 


••) T ; 


u = 


1 

2 


Pn 1 


df3 n ' 




Vn = 


1 „-l 9Hk 

2 rn Qa n i 


neZ 



L has the infinite matrix representation: 
/ n — > — oo 



ft 



n — > — oo 



o 



o 



-«n-l U n+ i 



n — > +oo 



n — > — oo 



o 



o 



-Vn-i o -w„+i 



V 



n — > +oo 



o 



o 



n — >■ +oo 



n — > — oo 



o 



o 



■ltn-1 U n+ i 



where w n = T r p n , v n = I> n , ft = f 
£2 norm: 



fci pi 
k 2 P2 



n — > +00 / 
We define the enstrophy norm which is the 



(a,/?) 5 



(IV.2) 
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Lemma IV. 1 The linear operator C maps £ 2 x £ 2 into £2 x ^2-' 



£ : l 2 xf 2 Hf 2 xf 2 



Proof: Notice that 



Pn — > \P\ , as \n\ — > 00 . 



Let 



= max \p n \ 



be the maximum of \p n \; then < 00, and from (IV. 1), we have 



|£n| < c I |a n +i| + \a n -i\ + \/3 n+1 \ + j , 



7n| < C j |d! n+ i| + |a n _i| + \Pn+l\ + |Ai-l| } , 



where c = \p* |T| k\ p 2 - h pi 



. Thus 



< 8c 2 



a 
P 



which proves the lemma. □ 

A complex number A is an eigenvalue of £, if there exists (a,j3) T 6 f 2 x f 2 such that 



(IV.3) 



Lemma IV. 2 The eigenvalues of the linear operator £ have the following properties: 

1. If A G C is an eigenvalue of C, then both A (the complex conjugate of X) and —X are 
also eigenvalues. 
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2. If X is a real eigenvalue, then X is a multiple eigenvalue. 

3. If X is a simple eigenvalue which is not real, then its corresponding eigenvector satisfies 
the relation (3 = ± ia. 

Proof: Since £ is a real linear operator, if A is an eigenvalue of C, then A is also an 
eigenvalue. Next we show that if A is an eigenvalue of C, then —A is also an eigenvalue. Let 
A and cu be an eigenvalue and a corresponding eigenvector, starting from the equation (II. 3), 
we have 



Xu k+n P = MP, k + (n- l)p) T uj k+(n _ 1)p 

+ A(-p, k + (n+ l)p) f u k+{n+1)p . 



Then, 
satisfies 



i-^Unp = A & k + (n- l)p) T u k+{n _ 1)p 

+ A(-p, k + (n+ l)p) f % +(n+1)p • 



Therefore, —A is also an eigenvalue. To prove claims 2 and 3, notice that the Hamiltonian 
7i k (II. 6) is invariant under the transformation 
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Therefore, if A is an eigenvalue and 




solves the linear system (II. 6), then 




also solves the system (II. 6). If A is real, then (a, f3) T can be chosen to be real. Assume that 

then rf = —1 which leads to a contradiction. Thus (a,(3) T and (— (3,a) T are linearly inde- 
pendent, and A is a multiple eigenvalue which proves claim 2. If A is simple and not real, 
then (a,j3) T and (— (3, a) T are linearly dependent and (IV.4) implies that (3 — ± ia, which 
proves claim 3. □ 

In fact, we have the following theorem. 

23 



(IV.4) 



Theorem IV. 1 The eigenvalues of the linear operator £ defined in (IV. 1), have the follow- 
ing properties: 

1. The eigenvalues of C are of four types: real pairs (c,—c), purely imaginary pairs 
(id,— id), quadruples (±c±id), and zero eigenvalues. 

2. If X is a real eigenvalue, then X is a multiple eigenvalue. If X is a simple non-real 
eigenvalue, then its corresponding eigenvector satisfies the relation (5 — ± ia. 

3- If J2k ^Apl = ®> then all the eigenvalues of C are either purely imaginary and in 
complex conjugate pairs (X = ic, —ic; c is real and not zero) or zeros. 

4- If \k\ = \p\, then zero is a multiple eigenvalue of C. 

5. If \k\ = \p\ and k + p D\ p \ (or k — p D\ p \), then all the eigenvalues for the system 
(III. 4) (or (III. 5)) are either purely imaginary and in complex conjugate pairs or zeros. 

Proof: Claims 1 and 2 are proved in Lemma IV. 2. Next we prove claim 3. If ^D\ P \ = 0, 
then If. is negative definite for (a,(3) T ef 2 x £2- In fact, there exists a constant 5 > 0, such 
that 



5 < -p n < 2 , Vn G Z . 



(IV.5) 



Thus, 




(IV.6) 
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Let A be an eigenvalue of C, and (a,(3) T be its corresponding eigenvector, which are written 
in terms of real and imaginary parts, 



\ = \ r + i \i, (a, [3) T = (a (1 \pV) T + i (a< 2 >, 0< 2 >) r . (IV. 7) 



Then both the real and the imaginary parts of e xt (a, [3) T , denoted by (a^\ (3^) and (c^ 2 - 1 , Z^ 2 - 1 )-! 
are real solutions to system (II. 6), where 



cr 



{1) = e Xrt \a w cos - a {2) sin A^l , 



a {2) = e Kt 



sin \d + cos \rf 



p{2) = e Mg(i) sinAji + ^) cosAti 



Denote by if (j = 1,2) the values of I~ k evaluated at (a (j \ (3 {j) ) T ] then 



I%\t) + I?\t) = e^Eft[(« (1) f + (« (2, ) 2 + (1) ) 2 + (^) 2 ] 
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4 1} (o) + / 2) (o) 



= EPn[(a {1) ) 2 + (a (2) ) 2 + (^ (1) ) 2 + (^ {2) ) 2 ] 



+ 0, 



by Lemma III. 1 and relations (IV.5, IV.6). Thus, 



e ZArt = 1 , i.e. A r = . 

This proves claim 3. To prove claim 4, notice that if \k\ = \p\, then system (II. 3) decomposes 
into systems (III. 4, III. 5, III. 6). Thus the vectors (a,P) T defined as 

a = 1 , a n = {n ^ 0) , (3 n = (VneZ); 



and 



A) = l, /^n = (ra^O), a„ = (VnGZ); 

give two linearly independent eigenvectors of C with eigenvalue zero. This proves claim 4. 
Claim 5 follows from the proof for claim 3 when restricted to system (III.4) or (III. 5). □ 
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Remark IV. 1 For a finite dimensional linear Hamiltonian system, it is well-known that the 
eigenvalues are of four types: real pairs (c, —c), purely imaginary pairs (id, —id), quadruples 
(±c±id), and zero eigenvalues [12] [13] [14]- There is also a complete theorem on the normal 
forms of such Hamiltonians [14]. 
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V. The Point Spectrum of the Linearized Two-Dimensional| 
Euler Equation: A Continued Fraction Study 



Rewrite the equation (II. 3) as follows, 

Pn & ^n— 1 

where ~z n = p n e in{e+ * /2) u~ k+np , 9 + 7 = tt/2, T = \F\e*, a = ±|r| 
\p\~ 2 - Let z n = e xt z n , where A G C; then z n satisfies 

where a n = A(ap„) _1 . Let w n = z n /z n -i [15]; then w n satisfies 

1 

a„ H = w n+1 . 

W n 

Iteration of (V.3) leads to the continued fraction solution [15], 



w n l) = a n _i + 



Pi h 
P2 h 



, Pn 



(V.l) 
\k + np\~ 2 — 

(V.2) 



(V.3) 



O-n-2 + 



(V.4) 



Qn-3- 



Rewrite (V.3) as follows, 



w n = 



—a n + w n+ i 

Iteration of (V.5) leads to the continued fraction solution [15], 



(V.5) 



w 



(2) 



a n + 



(V.6) 



a n +i + - 



n+2^ 



Before we study the two continued fraction solutions (V.4) and (V.6), we like to quote two 
theorems on the convergence of continued fractions [21]. 

Theorem V.l (Sleszyriski-Pringsheim's Theorem) The continued fraction 



K(a n /b n ) 



U l Q2 

°1 + h2 + ^^ 
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oo 

OO . 
n=l 



where {a n } and {b n } are complex numbers and all a n ^ 0, converges if for all n 

\b„\ > \a>n\ + 1 • 

Under the same condition 

|K(a n /6 n )| < 1 . 

Theorem V.2 (Van Vleck's Theorem) Let < e < ix/2, and let b n satisfy 

-tt/2 + e < arg{6„} < vr/2 - e 
for all n. Then the continued fraction K(l/6 n ) converges if and only if 

oo 

Notice that as n — > ±oo, 

a n — > a = — Aa" 1 |p| 2 . (V.7) 

Then we have the corollary. 

Corollary 1 If Re{a} ^ 0, or Re{a} = (]a| > 2), then the two continued fractions (V.4) 
and (V.6) converge. 

Proof: If Re{a} ^ 0, then there exists an positive integer TV and a positive constant e, such 
that 

—tt/2 + e < arg{a n } < n/2 — e 

for all |n| > N, or 

— n/2 + e < arg{— a n } < tt/2 — e 
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for all \n\ > N. In either case, applying Van Vleck's theorem, we have the convergence of 
the two continued fractions (V.4) and (V.6). If Re{a} = (\a\ > 2), then there exists an 
positive integer N such that 

\a n \ > 2 



for all |n| > N. Then applying Sleszyhski-Pringsheim's theorem, we have the convergence 
of the two continued fractions (V.4) and (V.6). □ 

Remark V.l In fact, as proved in Theorem VI. 5, Re{a} = and \a\ < 2 correspond to the 
continuous spectrum (= essential spectrum) of the system. 

When Re{a} ^ 0, or Re{a} = (\a\ > 2), as n — > — oo, 

ffi' 1 ' -> w {1) =a + = a + K(l/a) , (V.8) 

«;?> _> W W = 1 = _ K (l/a) . (V.9) 

° "r -i 1 ' 
a+— 

a + 

Both uA 1 ) and satisfy the equation, 

w 2 -hw- 1=0. (V.10) 
• When Re{a} ^ 0, the solutions of (V.10) can be written as: 



as n — > +oo, 



w ± = - 



± SVa 2 + 4 



(v.n) 



where 5 = sign(Re{a}) sign(Re{Va 2 + 4}). (Note that if Re{a} ^ 0, then RejVa 2 + 4} ^| 
0.) 
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When Re{a} = (\a\ > 2), let a = i£, £ is a real number, the solutions of (V.10) can 
be written as: 

w ± = l -k±5Je-4\ , (V.12) 



where 5 = sign(£) sign( v / ^ 2 — 4). 



Lemma V.l The solutions (V.ll) and (V.12) satisfy the inequality 



\w + \ > 1 > \w-\, 



and the continued fractions (V.8) and (V.9) have the values 



W W — w+ ^ ^,(2) — w _ 



Proof: First we show that > 1 when Re{a} ^ 0. 

w + w^ = ^(ja| 2 + |a 2 + 4| + dl5 \J~a 2 + 4 + aS^a 2 + . (V.13) 

Let a = ai + ia 2 , \/a 2 + 4 — bi + ib 2 \ then 

hb 2 = ai a 2 , (V.14) 

b\-b\ = a\-a 2 2 + A. (V.15) 

From (V.14), we have 

a\{a 2 b 2 ) = bKa^i) . 
Thus, a 2 b 2 is either zero or of the same sign as a\b\. Therefore, 



~a^~o? + 4 + a5\/a + 4 = 2<5[aA + a 2 b 2 ] > . (V.16) 
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Together with 

\a 2 + 4| > 4- \a\ 2 , 

we have \w+\ > 1. When Re{a} = (|a| > 2), it is obvious that > 1. Notice that 
w + w_ = — 1, we have |u>_| < 1 in both cases. Thus, we have 

\w + \ > 1 > . 

From the relation w + w^ — — 1, when Re{a} ^ 0, Re{w + } and Re{tu_} are of opposite signs. 
When Re{a} ^ 0, Re{u>W} and Re{w^} are of opposite signs, and Re-jW 1 )} and Re{w + } 
are of the same sign; thus, = w + and = u>_. When Re{a} = (|a| > 2), by 
Sleszyriski-Pringsheim's theorem, 

|K(l/o)| < 1 ; 

then, 

= \a + K(l/a)| > \a\ - |K(l/a)| > 1 . 
Thus, uA 1 ) = u> + and —W-.D 

Definition 5 Define as follows 

= for n < 1, = for n > 1. 

Then ujM solves (V.3), provided that = i.e. 

f = a 0+ ( ■ 1 - ) + ( 1 - ) = , (V.17) 

a_ 3 + z 03 + 

where / = f(\, k,p), A = A/a. Let satisfy 

„..(*) _ _(*) /_(*) . 
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then asn-> +00, 
and asn-> —00, 

4* } ~ (^+) n • 

Thus by Lemma V.l, z^ G £ 2 . Therefore equation (V.17) determines eigenvalues. 

If fl D\ p \ = 0, then p n < for any n G Z. If Re{A} 7^ 0, then Re{a n } 7^ and are of a 
fixed sign for any n E Z. Then Re{/} 7^ 0. Therefore, in such cases, there is no eigenvalue 
with nonzero real part. This fact is already obtained in Theorem IV. 1. Moreover, we have 
the following fact. 

Lemma V.2 //S^ fl D\ p \ = 0, then equation (V.17) determines no eigenvalue. 

Proof: As discussed above, if fl D\ p \ = 0, then the possible solution A to (V.17) has to be 
imaginary. Therefore, a has to be imaginary. Rewrite equation (V.2) as follows 

Lz n = —\z n+ i - z n _i] = az n . (V.18) 
P 

If n/)|p| = 0, then < p n j ' p < 1 for all n G Z. Thus, ||L|| < 2. Then if A is an eigenvalue, 
|S| < ||L|| < 2. Notice that equation (V.17) should be solved under the condition Re{a} 7^ 
or Re{a} = (\a\ > 2); thus, in this case, equation (V.17) determines no eigenvalue. □ 

Remark V.2 In fact, as proved in Theorem VI. 5, Re{a} = and \a\ < 2 correspond to 
the continuous spectrum (= essential spectrum) of the system. Thus, fl D\ p \ = 0, the 

point spectrum is empty (cf: Theorem VI. 5). Then the problem is reduced to solving equation 
(V.17) under the conditions fl D\ p \ 7^ 0, Re{a} 7^ or Re{a} = (\a\ > 2). 

Example: Let p = (1, 1) T , in this case, only one class labeled by k = (1,0) T has 
no empty intersection with D\ p \ (the other class labeled by k = (0, 1) T gives the complex 
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conjugate of the system led by the class labeled by k = (1,0) T ). For this class, \p n j p\ < 1 
for all n E Z. Thus, the linear operator L defined in (V.18) has norm ||L|| < 2. Therefore, 
equation (V.17) determines no real eigenvalue. Numerical calculation on equation (V.17) 
gives the eigenvalue: 

A = 0.24822302478255 + i 0.35172076526520 . 

By Theorem IV. 1, equation (V.17) determines a quadruple of eigenvalues, see figure 3 for an 
illustration. 
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VI. The Spectra of the Linearized Two- Dimensional 
Euler Equation: An Infinite Matrix Study 

VI. 1. The General Setup 

Rewrite (II. 3) as follows: 



z n = la 



Pn-l z n~l + Pn+1 z n+l 



, (neZ) 



(vi.i; 



where 



r = |r|e i7 , # + 7 = tt/2, a = - |r| 



p 1 kj 

P2 h 



Relabel {z n } as follows: 



z n — z 2n i n > 1 , 

Z-n — z 2n+l , Tl > 0] 



then 



z 2n = ia 



Pn-l z 2(n-l) + Pn+1 z 2(n+l) 



(n > 2) 



(VI.2) 



z 2n+l — %a 



P-n+l ^2(n-l)+l + P-n-1 z 2(n+l)+l 



, (n>l) (VI.3) 
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Z2 



ia 



Po zi + p 2 z 4 



Z\ = ia 



Pi z 2 + p-i z 3 



(VIA) 



(VI.5) 



for z = (zi,z 2 ,- ■ -) T - Notice that Equations (VI. 2) and (VI. 3) are decoupled, the coupling 
between components of z with even and odd indices is through Equations (VI. 4) and (VI.5). 
The right hand side of (VI.2)-(VI.5) define a bounded linear operator Ca ■ (2 ^ ^2, with 
the infinite matrix representation, 



ia 



( 





Pi 


p-1 














\ 




Po 








P2 















Po 











P-2 
















pi 











P3 















p-1 











P-3 




\ 





















(VI.6) 



More importantly, 



P = - \P\ 



1-2 



as ra 



00 . 



(VI.7) 



Define the infinite matrix 



B = ib 



( 1 1 

1 1 

1 1 O 

1 1 

1 1 



(VI.8) 



V O 



o 
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where b = ap = —a \p\ . Define the infinite matrix C as 



C = A-B, 



(VI.9) 



that is, 



C = ia 



( ° 

Po 

Po 





Pi 



pi 





p-i 
p 2 



p 3 o 
P-i o p_ 3 



o o p_ 2 o o o 



\ 



(VI.10) 



V 



o 



o 



where p n = p n — p. Denote by Cb and Cc the bounded linear operators with the infinite 
matrix representations by B and C. According to Duren [18], La, Lb and Cc are called 
(2x2 + l)-operators, since their entries c„ in+m satisfy the condition c„ in+m = if \m\ > 2. 
Cb is a (2 x 2 + l)-operator with constant coefficients, since its entries c„ in+m is independent 
of n when n > 2; and is self-adjoint. 

Theorem VI. 1 The bounded linear operator Cc : £2 ^ ^2 is a compact operator. 

Proof: Denote by C^ the linear operator represented through the matrix Cn x n ob- 
tained from C by replacing its entries c„ hn by 0, when m > N. Let I | be a bounded 




sequence in £ 2 ; then {C^z^ 



is a bounded sequence in which each element has only one 



nonzero component, i.e. 



(C i c ) z U) ) n = 0, whenn>l. 
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Thus, there exists a subsequence {z^}, such that {C^z^} converges in £ 2 . Similarly, we 
can get a subsequence of {z^}, denoted as {z^}, such that {C^z^} converges in £ 2 , 
and so on. Therefore, we have a nested list of subsequences: 



z (ii) z (ia) 

z (2l) z (2 2 ) 



We choose the subsequence {z^ nn ^} of {z^}, which is the diagonal of the above list. There 
exist constants ( and N Q , such that 



c 



< St; , for all n > N and all n . 

n 2 



(VI.11) 



For any e > 0, choose N large enough, such that 



iV 2 3 



(VI. 12) 



Since the subsequence converges, there exists N, such that 



<ge, Vji,32>N. 



(VI.13) 



Let iVi = max {iV, iV"|; then 
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< - e Vn , n > Ni 



(VI. 14) 



Thus 



C c - Cc z^ 



a 



< 



+ 



C 1 



c 



+ 



C^z™ - C c z^ 



Therefore, ^Cc 2^™")} is a Cauchy sequence in £ 2 ; thus converges. This proves that Cc is a 
compact operator. □ 

Remark VI. 1 In fact, a theorem of Achieser and Glasmann [22] [18] states that a (2M+1)- 
operator is compact if and only if its diagonal sequence entries tend to zeros, i.e. c ni „ +m — > ; 
as n — > oo for each fixed m, \m\ < M. Here we give the proof for self-containedness. 

VI. 2. The Spectra of the Linear Operator Cb 

Next we will follow a theory of Duren [18] to study the spectra of the constant-coefficient 
infinite-matrix bounded self-adjoint operator 
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The characteristic polynomial for the difference equation 



(B-\I)z = 0, (VI.15) 



where I is the identity matrix, is defined as: 



Jb(w, X) — ib - Xw +ib w . (VI. 16) 



Define the rescaled characteristic polynomial as follows: 



fs(w, X) — 1 — Xw +w\ (VI.17) 



where A = ibX. In fact, /b(u>, A) is the characteristic polynomial for the difference equation 



(B-XI)z = 0, (VI. 18) 



where B = —ib 1 B. The roots of A) are: 



tu* , -tu* , — , — , (VI.19) 



where 
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1/2 



A + \ \ 2 - 4 



Definition 6 The spectral curve of the linear operator (with the infinite matrix rep- 
resentation by B), denoted by Cq, is defined to be the set of all A e C for which the char- 
acteristic polynomial f B (w,X) has a root of modulus one. The spectral point-set of the 
operator C B , denoted by P B , is defined to be the set of all A e C for which the characteristic 
polynomial f B (w,X) has a multiple root. Denote by S B (\) the number of roots of f B (w,X), 
of modulus less than 1 (counted with multiplicity). 

Notice that 

A = vol + ™*~ 2 . (VI.20) 

Let if* be a root of modulus 1 (then all the four roots are of modulus 1), u>* = e' t9 , 9 e [0, 27r); 
thus the spectral curve C B is the segment of the real axis, 

C B : A = 2cos2#, 9e[0,2n). (VI.21) 



See Fig. 4. The spectral point-set P B consists of two points, 

P B : A = ± 2 , (VI.22) 
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which are the boundary points of the spectral curve C B . At A = ± 2 , the four roots of 
f B (w, A) are, 



at A = 2 : 1,-1,1,-1; 
at A = — 2 : % , —i , —i , % . 



The function S^(\) is 



, if A e C § , 

S 5 (A) = < ' (VI.23) 

2, if \(£C § . 



The general solution to (VI. 18) is 



z n = c x v% + c 2 {-w*) n + c 3 w~ n + c 4 (-w*)- n , (VI.24) 
if A£P B , 

z n = ci«j" + c 2 n + c 3 (-M) t ) n + c 4 n (-«;»)" , (VI.25) 
if A G (then w* — 1 , i) ; 



under the restrictions: 



-A *i + z 2 + * 3 = , (VI.26) 
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z i — X z 2 + £4 = . 



(VI.27) 



Lemma VI. 1 The general solution (VI. 24, VI. 25) is in £ 2 if and only if \w*\ < 1 and 
c 3 = c 4 = in (VI. 24) or \w*\ > 1 and c\ = c 2 = in (VI. 24). 

Proof: See ([18], pp. 24, Lemma 5). □ 

Definition 7 Let C : £ 2 l— ► ^2 ^ e fl linear operator. The set of points cr p (jC) in the complex 
X-plane C such that (C — XI) has no inverse (i.e. £ — XI is not 1-1), is called the point 
spectrum of C. The set of points a r (£) in C such that (£ — A/) -1 exists and is a linear 
operator with domain not everywhere dense is called the residual spectrum of C. The set 
of points cr c (C) in C such that (C — XI)' 1 exists and is an unbounded linear operator with 
domain everywhere dense is called the continuous spectrum of C The set of points p{£) in 
C such that (£ — A/)" 1 exists and is a bounded linear operator with domain everywhere dense 
is called the resolvent set of C. The set a(C) = cr p (jC) U o>(£) U a c (C) is called the spectrum 



Without loss of generality, assume < 1. Then A is an eigenvalue if and only if there 
exists a non-trivial solution (ci,c 2 ) to the following system, 



ofC. 




I -Aw* + w\ + w\ Aw* + w\ - vol \ 



(VI.28) 
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Theorem VI. 2 The £ 2 point spectrum cr p (B) of the linear operator C B is empty. 



Proof: The determinant 



det 



/ -Au>* + vol + wl Au>* + vol - vol \ 



\ - \vol + vol -to* - \vol + wf J 



2voi 



wl-2\vol + (\ 2 -l) =0, 



(VI.29) 



implies that 



vol - 2Xwl + A - 1 = , (VI.30) 



since vo* ^ 0. vo* is a root of /b(u>, A) (VI. 17), 



<-A< + l = 0. (VI.31) 



From (VI.30, VI.31), we have 



vol = ^j 1 ■ (VL32) 



Notice also that 



7 o q A^ — 2 A 
A = vol + w„ = — - — + 



A A 2 - 2 ' 



44 



which implies that A = ± 2 . Then \w*\ = I. Thus if \w*\ < 1, then Equation (VI. 28) has 
only trivial solution. Therefore, the point spectrum of C B is empty; equivalently, the point 
spectrum of C B is empty. □ 

Theorem VI. 3 The £2 residual spectrum <J r (B) of the linear operator C B is empty. 

Proof: A G cr r (B) if and only if the dimension of the orthocomplement of (C B — XI) o £ 2 
is positive and (C B — A/)" 1 exists. From the inner product relation 

since C B is self-adjoint, C B = C*^ (the adjoint of £§), where ( , ) denotes the inner product 
over the complex field, we have that if A G cr r (B), then A G o- p (B). By Theorem VI. 2, c p (B) 
is empty; thus, a r {B) is empty; equivalently, a r (B) is empty. □ 

Since iC B is self-adjoint, this theorem is well-known, but we furnish a short proof here. 
From (VI. 15, VI. 18) and (VI. 21, VI. 22), the spectral curve C B for the linear operator C B is 
the segment of the imaginary axis, 

C B : A = i2b cos 20, 6 G [0, 2ir) ; (VI.33) 

the spectral point-set P B for the linear operator C B is 
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P B : A = ± i2b , 



(VI.34) 



which are the boundary points of the spectral curve Cb- 

Theorem VI. 4 The £ 2 continuous spectrum a c (B) of the linear operator Cb is the spectral 
curve, cr c (B) = Cb ■ The £2 resolvent set p(B) of the linear operator Cb is the complement 
of Cb in the finite complex plane C , p(B) = (Cb)' ■ 

Proof: First we show that if A G C^, then A G cr c (B). Since both o~ p (B) and o- r (B) are 
empty by Theorems VI. 2 and VI. 3, for any A G C B , (£g — A/) -1 exists and is everywhere 
densely defined. We need to show that (C^ — A/) -1 is unbounded. For any A G C#, there 
exists a root of /b(u>, A) of modulus one, 

w, = e i6 , Be [0, 2tt) . 

Define the elements 



'ind 



n<N, 



0, n > N. 



Then e £ 2 for each finite N, and || z^\\ — > 00, as N — > 00. There exists a constant d 
independent of N, such that 



(B-XI)z^ <d, ViV 
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Thus 





Z (N) 






-XI) 





oo , as N — > oo . 



Therefore, (£g — XI) 1 is unbounded, and A G (T C (B). Next we show that if A ^ C^, then 
A G p(-B). For any A ^ C^, the corresponding roots of /b(u>, A) are (VI. 19), such that 



luU = I— wJ < 1 < 



(VI.35) 



For any i/ G f2, we want to construct a solution to 



(VI.36) 



using the method of variation of coefficients. Explicitly, we need to solve 



z n — X z n+ 2 + z n+ 4 — y n +2 , {n > 1) 



(VI.37) 



under the constraints 



-Xzi + z 2 + z 3 = y 1 , 
k - Xz 2 + Zi = y 2 . 



(VI.38) 
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Assume a solution to (VI. 37) has the form 

z n = ci 1} < + c%\-w*) n + c^w; n + c£\-w,)- n . (VI.39) 

If 

Ac« < +1 + Ac^ (-w*) n+1 

+Acf + Aci 4) (-«;,) - {n+1) = , (VI.40) 

Ac« <+ 2 + Aci 2 ) (-^)« +2 

+Acf W ;(" +2 ) + Aci 4) (-«;,)- ( " +2) = , (VI.41) 
Ac« <+ 3 + Ac( 2 ) (-w*) n+3 

+Acf ^;( n+3 ) + Aci 4) (-w*)~ (n+3) = , (VI.42) 
AcW <+ 4 + Ac?* (-wX +4 
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n+2 j 



(VI.43) 



where Acf = - c%\ (£ = 1,2,3,4), then z n given in (VI.39) solves (VI.37). Solving 
(VI.40- VI.43), we have 



where 



n(<) 

AcW = (-l)Sn +2 -^, (€=1,2,3,4) 



(VI.44) 



W„ 



w 



n+1 



W 



n+2 



{-w.) n+1 w: {n+1) (-^)-(" +1 ) 

{-w*) n+2 w7 {n+2) (-w*)~ (n+2) 



< t+3 (-w») n+3 w*~ (n+3) (-w,)-^ 3 ' 



-(ra+4) 



(VI.45) 



w*" (n+1) 



-w*)~ (n+1) 

(-w*) n+2 w ; in+2) (-^)~ (n+2) 

(-w*) n+3 uC (n+3) (-w*)~( n+3 ) 



2w; 



uC 4 - 1 



(VI.46) 
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w 



n+1 w: in+1) 



-(n+l) 



,n+2 „„-( n + 2 ) 



"+ 3 w ," (n+3) 



-(n+3) 



2(-w*) _n 



1 — wl 



(VI.47) 



n+l 



W 



w? +2 {-wS) n+2 



-(ra+2) 



n+3 (_^)«+3 (_^)-(n+3) 



(VI.48) 



M 4) 



»+2 (_^)-+2 ^+ 2 ) 
»+3 (_^)«+3 ^ (n+3) 



it; 



w 



2(-w, 



l-wt 



(VI.49) 



The W n defined in (VI. 45) satisfies the Wronskian relation 
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w n+1 = w n . 



(VI.50) 



The representation (VI. 44) can be extended to n > 0. Choose cjp = 0, we have 



n— 1 i n— 1 

cS? = E Acf = 77T E ^ (-1/ I/i + 2 , (n > 1) . (VI.51) 

i=0 W j=0 



The expressions (VI.46-VI.49) lead to 



c (2) = 



2[1 E^% + 2, (VI.52) 



Wo j=0 



^ E(-^)" j % + 2, (VI-53) 



W 



3=0 



c (3) 



ci 4 ) = 



r ' E^% + 2, (VI.54) 



21 "'' ] £< >>■■)' <Jj-2- (VL55) 



, =0 



With these representations of z n given by (VI. 39) is a special solution to (VI. 37). The 
general solution to (VI. 37) is 



z n = (<£> + aW) < + (ci 2 ) + a< 2 >) (-«,,)' 
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+ (cf + a^) w: n + (cj? + a®) (-w*y n , (VI.56) 



where (£ = 1,2,3,4) are arbitrary constants. Set 



a® = 



2 1 ~ wj] ~ 

777 2^wi Vj+2, 

Wo j=0 



(VI.57) 



,(4) 



2 [!-«*] ^, v - 



i=o 



(VI.58) 



Let 



U = + ci 2 )(-^r + (cf + a (3) K~" + (#> + a( 4 ))(-^)-" ; (VI.59) 



then from expressions (VI.52-VI.55) and (VI. 57- VI. 59), we have 



fn 



2 [1 - w, 
W 



-41 n-1 



x: [<- j + (-™*r 



%+2 



i=o 



%+2 



j=n 



(VI.60) 



Finally, 



z n = aWw: + aW{-w.) n + f n . 



(VI.61) 



52 



Next we choose and to satisfy the constraints (VI. 38): 



M 



,.(2) 



Vi + A/i -f 2 -fa 
V2 - fi + A/ 2 - / 4 



where 



M = 



-Aw* + w 2 + wf Aw* + w 2 — wf 



w* - Aw 2 + wf -w* - Aw 2 + w * 



As shown in the proof of Theorem VI. 2, M is nonsingular. Then, 



,.(2) 



M" 1 



yi + A/i - / 2 - h 
V2 - /i + A/ 2 - / 4 



where 



M~ l = K 



-w* 



Aw 2 + w? —Aw* — w 2 + w; 



— w* + Aw 2 — wf —Aw* + vol + wf 



where « = 2w^ (w* - 2 Aw 2 + A 2 - l) . Thus 



w* , (-w* 



M -1 



l/i + A/i - h~ h 

V2-fl + A/2 - /4 



solves Equation (VI. 36). Rewrite / n given in (VI. 60) as follows: 
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fn = J2d(n,j) Vj, 

3=1 



(VI.65) 



where 



9(n,j) = < 



0, 



mgfl [w n- j+ 2 + (_^)n-i+2] , 2 < j < n + 1 ; 



^ kr"- 2 + (-^) j -"- 2 l , i > n + 2 . 



Rewrite z n given in (VI. 64) as follows: 



3=1 



(VI.66) 



where 



G(n,j) 



5 1 , j + ~Xg(l,j)-g(2,j)-g(3,j) 
5 2d -g(l,j) + \g(2,j)-g(4,j) 



+ g(n,j) , 



(VI.67) 



where Sgj is the Kronecker delta: Sgj = 1 (i — j), Sej = (i ^ j)- From the expression 
(VI.67), we see that there exists a constant X independent of n, j; such that 



J2\G(n,j)\ < K , V n = l,2, 
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(VI.68) 



oo 

J2\G(n,j)\ < K, Vj = l,2,- 

n=l 



(VI.69) 



Then, we have the norm relation, 



klloo = sup \z n \ < sup ^2\G(n,j)\ \yj\ 

n j=1 



< 



sup ^2\G(nJ) 



(VI.70) 



and the l\ norm relation, 



blh = lim 

iV^oo 



" N 




E \ Z n\ 


< lim 




iV— >oo 


_n=l 






= lim 




JV->oo 



N oo 
n=l jr' = l 



oo JV 

El%'lElG(n,j)l 

j=l n =\ 



< K J2\y J \=K\\y\\ 1 . 

i=i 



(VI.71) 



Thus the linear operator defined in (VI. 66) which maps y into z, is bounded in and 
l\. Therefore, by Riesz convexity theorem [23] [24] [18], (£g — A/) -1 defined in (VI. 66) is 
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bounded in £ 2 . Since by Theorems VI. 2 and VI. 3, (L B — XI)^ 1 exists and is everywhere 
densely defined and is bounded, we have A G p{B). In summary, we have shown that if 
A G C B , then A G a c (B); and if A ^ C B , then A G p(-B); thus, cr c (S) = C B and p(5) = (C^)'. 
Equivalents, a c (B) = C B and p(B) = (C B )'. □ 

In summary, the spectrum of Lb is as depicted in Figure 5. 

VI. 3. The Spectra of the Linear Operator La 

Now we apply Weyl's essential spectrum theorem [19] to obtain the spectral theorem for La- 

Theorem VI. 5 (The Spectral Theorem of La) 1- U^k n ApI = ®> ^ en ^ e en ^ re 
£ 2 spectrum of the linear operator La is its continuous spectrum which is the spectral 
curve Cb defined in (VI. 33), i.e. ct{La) = <t c (Ca) = Cb- See Figure 5. 

2. If fl D\ p \ 7^ 0, then the entire essential £ 2 spectrum of the linear operator La is its 
continuous spectrum which is the spectral curve Cb defined in (VI. 33), i.e. c ess (LA) = 
Cc(La) = Cb- That is, the residual spectrum of La is empty, a r (LA) = 0- The point 
spectrum of La is symmetric with respect to both real and imaginary axes. See Figure 
6. 

Proof: First, we want to show that in both cases, the residual spectrum of La is empty. 
By Weyl's essential spectrum theorem [19], the essential spectrum of %La is the same with 
the essential spectrum of iLb, cr ess (i£ J 4) = cr ess (i£ B ) = iC B . Let i\ r G o- r (iLA), then 
i\ r G iC B . By the argument in the proof of Theorem VI. 3, i\ r G cr p ((iCA)*), where {iLa)* 
is the adjoint of iLa, 

{iLaY = iL B + (iLc)* ■ 
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By Weyl's essential spectrum theorem [19], the essential spectrum of (iCa)* is the same with 
the essential spectrum of iCs, <?ess((i£A)*) = cress {i£>Ei) — iCb- Thus, i\ r G cress a)*)- 
Since cr ess ((i£ A )*) and (t p {{iCa)*) are disjoint, a r {iCA) = 0- The claim <t p (£a) = in case 
1 follows fom the proof of Lemma V.2 and the fact that the spectral curve C B corresponds 
to Re{a} = and \a\ < 2. The property of o- p (£a) in case 2 has been proved in Theorem 
IV. 1. Then Weyl's essential spectrum theorem implies the rest of the claims. □ 

Remark VI. 2 By the above theorem, the computation of eigenvalues is reduced to the case 
that E^,nZ)|p| 7^ 0. By Corollary 1, if A ^ C# = o- eS s(C-A), the two continued fractions (V.4) 
and (V.6) converge. And solutions of equation (V.17) lead to eigenvalues. 

Remark VI. 3 The width of the continuous spectrum <t c (£a) is A\b\, where b = — a\p\~ 2 and 
a = \ |r| ^ 1 - 1 . Although \a\ can increase to infinity as \k\ increases to infinity, a is 
essentially a scaling factor for £ a (is can be seen in the expression for the infinite-matrix A. 



Next we discuss an alternative way of representing eigenvalues. This approach is not useful 
for practical computation. Consider the linear difference equation, 

(A-XI)z = 0, (VI.72) 
where A defined in (VI. 6) is the representation matrix of La- Explicitly, 



Pn-l Z 2 ( n -1) - A Z 2n + Pn+1 ^2(n+l) = , (n > 2), 

. P-n+1 z 2(n-l)+l — ^ Z 2n +1 + P-n-1 z 2(n+l)+l = , (n > 1), 



(VI.73) 



under the constraints 
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-A z 1 + pi ^ 2 + p-i z 3 = , 

< 

k p 21 - A z 2 + p 2 z 4 = , 



(VI. 74) 



where A = {ia)~ l A. By Theorem VI. 1, the linear operator La is a compact perturbation of 
L B - Thus, the difference equation (VI. 73) is of Poincare-Perron type. The Poincare-Perron 
theorem stated specifically for the difference equation (VI. 72) is as follows [16], [17], [18]: 

Theorem VI. 6 (Poincare-Perron Theorem) For any A G C , let w*, — u>* ; — ^j- 
be the roots of the characteristic polynomial Jb(w, A) defined in (VI. 16), which are given in 
(VI. 19), where \wj\ < 1. Then there exists a fundamental set of solutions z$ (j = 1,2,3,4) 
to the difference equation (VI. 73), such that 

= (.7 = 1,2), 

(VI.75) 

= Ar, (J=3,4). 
\w*\ 

It is easy to see that 

z®et 2 , (j = 1,2); z®(£l 2 , (j=3,4); 

if |u>*| < I. By definition, when \ C B (defined in (VI. 33)), |u>*| < 1. Next we study the 
conditions for the point spectrum of La- Let 
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lim sup 



xCi) 



lim sup 



xCi) 



(VI.76) 



Substitute z n into the constraints (VI. 74), we have 



(VI.77) 



where 



M 



— AZ 1 + p\Z 2 + P-lZ 3 —AZ 1 + piZ 2 + P-lZ 3 
\ P0?i — AZ 2 + P2Z4 PqZ 1 — AZ 2 + P2^4 / 



(VI.78) 



Theorem VI. 7 If A ^ Cb (the spectral curve for Cb, defined in (VI. 33)), and detM = 
(where M is defined in (VI.78)), then A G cr p (A) (the point spectrum of La)- 

Proof: If det M — 0, then there is a nontrivial solution to (VI.77). Thus there is a nonzero 
solution to (VI. 73), which satisfies the constraints (VI. 74). Therefore, A is an eigenvalue. □ 
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VII. Conclusion 

In this paper, we study the linearized two-dimensional Euler equation at a stationary state. 
This equation decouples into infinite many invariant subsystems. Each invariant subsystem is 
shown to be a linear Hamiltonian system of infinite dimensions. Another important invariant 
besides the Hamiltonian for each invariant subsystem is found, and is utilized to prove an 
"unstable disk theorem" through a simple Energy-Casimir argument. The eigenvalues of 
the linear Hamiltonian system are of four types: real pairs (c, — c), purely imaginary pairs 
(id,— id), quadruples (±c±id), and zero eigenvalues. The eigenvalues are studied through 
continued fractions. The spectral equation for each invariant subsystem is a Poincare-type 
difference equation, i.e. it can be represented as the spectral equation of an infinite matrix 
operator, and the infinite matrix operator is a sum of a constant-coefficient infinite matrix 
operator and a compact infinite matrix operator. We have a complete spectral theory. The 
essential spectrum of each invariant subsystem is a bounded band of continuous spectrum. 
The point spectrum can be computed through continued fractions. 

This study is the first step toward understanding the unstable manifold structures of 
stationary states of the two-dimensional Euler equation, which we believe to be the key for 
understanding two-dimensional turbulence. In particular, we will be interested in investigat- 
ing whether or not the unstable manifolds of 2D Euler equations are degenerate (i.e. figure 
eight structures). Degeneracy will imply that the dynamics of 2D Euler equations is not 
turbulent. 
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K{A} 



Figure 3: The quadruple of eigenvalues determined by equation (V.17) for the system led by 
the class labeled by k = (1, 0) T , when p = (1, 1) T . 
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Figure 4: The spectral curve . 
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Figure 5: The continuous spectrum of Lb and La- 
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Figure 6: The spectrum of Ca in case 2. 



Figure 1. Caption: An illustration on the locations of the modes (k' = rk) and (\k'\ = \k\) 
in the definitions of El and E\ (Proposition 1). 

Figure 2. Caption: An illustration of the classes and the disk D\ p \. 

Figure 3. Caption: The quadruple of eigenvalues determined by equation (V.17) for the 
system led by the class Ejg. labeled by k = (1, 0) T , when p = (1, 1) T . 

Figure 4. Caption: The spectral curve . 

Figure 5. Caption: The continuous spectrum of Cb and La- 

Figure 6. Caption: The spectrum of C A in case 2. 
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